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We use the Feynman-Kleinert variational method [R.P. Feynman and H. Kleinert, Phys. Rev. 
A 34, 5080, (1986)] to calculate the partition function and effective pinning energy of a single flux 
line near a columnar pin in type II superconductors. It is found that there is a phase transition 
between a low temperature phase where the flux line is localized near the columnar pin and where 
its internal modes fluctuations are bounded, and a high temperature, depinned phase where the flux 
line is delocalized and its internal fluctuations are those of a free line. 



I. INTRODUCTION 

During the past decade, the static and dynamic prop- 
erties of flux lines in high temperature superconductors 
(HTSC) in the presence of a random array of pinning 
centers have been the subject of a considerable amount 
of work (for recent reviews, see Blatter et alE, and Nat- 
termann and SchcidlEI) . Among the various types of pin- 
ning phenomena which have been investigated, a particu- 
larly promising one is the relatively strong pinning which 
results rfppm the interaction of flux lines with columnar 
defectsraa. The latter can be viewed as long, parallel 
damage tracks of diameter do ~ 50 — 70A, produced by 
heavy ion (e.g. Sn, Pb) irradiation at high (~ GeV) en- 
ergies. Such columnar defects were found experimentally 
to lead to a significant enhancement of the irreversibil- 
ity line, with critical currents increasing by as much as 
three orders of magnitude at T — 77 K in thallium- 
based compoundsB. On the theoretical side, much of 
our present understanding of the pinning of flux lines by 
these columnar defects originates in the wotk-of Nelson 
and Vinokum, who used the boson analogyutl to study 
the physics of individual flux lines as well as flux line 
lattices and liquids in the presence of columnar pins. 

In the present paper, we make an attempt at describ- 
ing the physics of a single flux line-columnar pin pair, 
without resorting to the boson formulation. The flux 
line will be parametrized by the three-dimensional vector 
R(z) = (r(z),z), where the two-dimensional vector r(z) 
describes the position of the flux line element at height 
z in the (x, y) plane. We assimilate the flux line to an 
elastic string whose internal degrees of freedom can be 
described by the HamiltonianQ 

r L I /dr\ 2 

H - = l *2*(e) <"> 

where L is the sample thickness in the z direction and E\ 
is the tilt energy per unit length of the flux line. For a su- 
perconductor with uniaxial symmetry, assuming that the 
external magnetic field H||c, e"\ can be approximated by 
~ £ 2 £okiK, where e = A a (,/A c is the ratio of the London 



penetration depths along the (ab) plane and the c axis, 
£o = (0o/47rA a h) 2 (<po = hc/2e is the flux quantum), and 
K = A a b /£ a ;, is the ratio of the London penetration depth 
to the coherence length in the (ab) plane. 

The total Hamiltonian of our system of one single flux 
line in the presence of a columnar pin is then given by 

H = L dz h~<Tz) +v{ ^ (L2) 

where the z-independent function V(r) denotes the pin- 
ning potential due to the columnar defect. V(r) will 
be taken as a very localized function, V(r) — Uq 
for |r| < b and V(r) « for |r| > bo, withO b = 
m&x(d /2,V2U) and U Q » e ln[l + (d /2 % /2^ Qb ) 2 ]/2. 
If we take as the origin of free energies the free energy 
of a free flux line (i.e. a flux line far from the columnar 
pin), then the partition function of our system can be 
written as the normalized path integralQ 

/ Vv(z) exp { - f3 J Q L dz [f (g) 2 + V(v(z))] } 
JVr(z) exp{-/3/ L dz^(g) 2 } 

(1.3) 

with P = 1/kgT the inverse temperature (kg is Boltz- 
mann's constant). In its simplest form, which happens 
also to be the form relevant to our present problem, the 
boson analogy consists in exploiting the formal corre- 
spondence between the above, normalized partition func- 
tion, and the partition function of a quantum particle of 
mass m in imaginary timaijli3(= r), as can be seen by 
using the replacementsa z — > r, kgT — > h, ix — > m and 
L — ► f3H (h is Planck's constant h divided by 2tt). This 
correspondence is then used to express the partition func- 
tion Z as 

Z = T J dr a J dr b p(r a ,r b ; L) (1.4) 

where T is a normalization constant with dimension 
(length)^ 2 , whose precise value is unimportant to the 
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arguments that follows, where r a = r(0) and = r(L) 
denote the endpoint positions of the flux line and where 
p(r a , Tb] L) is the density matrix 



(r a |e-"^ T |r 6 > 



p(r a ,r 6 ;L) 

associated with the "quantum" Hamiltonian 
(k B T) 2 



H = 



2ei 



V(v) 



(1.5) 



(1.6) 



(here Vj_ = (d x ,d y ) is the gradient in the (x,y) plane). 
Using a complete set {ip n (r)} of eigenf unctions of the op- 
erator TL with eigenvalues E n , the density matrix can be 
written in the form 



p{v a ,r b ;L) =J2 e ~ LEn/kBT ^ r ^n(n) 



(1.7) 



Equations (|l.4[ ) and (1.7) lead to the following expression 
for the partition function 



dr a ip n (r a ) / dr b ipKn) (1. 



It should be noted that, due to the appearance of temper- 
ature T in the quantum Hamiltonian TL, both the eigen- 
values E n and eigenfunctions tjj n of TL will depend on 
temperature in a rather complicated way. At low enough 
temperatures however, we expect TL to possess many low 
lying, localized bound states. Under this assumption, we 
see that in the limit L — > oo, the above sum is dominated 
by the term corresponding to the ground state E < of 
the quantum Hamiltonian TL, and may be approximated 
by 

z~r e - LEo/kBT j ' dv a M*a) fdn^in) (1.9) 

Taking — Eq as an approximate expression for the effec- 
tive pinning energy U e (T) per unit length of the flux line 
trapped in the,-cplumnar pin, simple quantum mechan- 
ical arguments^ lead to to the conclusion, that in this 
low-temperature regime, U e (T) is given byo : 



U a {T) = -Eo*U 



ci(k B T) 2 



(1.10) 



(here c\ is a numerical constant of order unity) . We thus 
see that thermal fluctuations of the flux line average out 
the pinning potential of the columnar pin, reducing its 
strength from the bare value Uq to U e (T) < Uq. Defining 
the characteristic temperature T* such that Eq(T*) ~ 0, 



k B r* ~ b Q ^Fyu , 



(1.11) 



it is found that two distinct physical regimes will emerge, 
depending on the value of the temperature with respect 
to the characteristic value T*, namely : 



(i) for T <C 7 1 * , U e (T) is reduced from its bare value ac- 
cording to (here u(z) = r(z) — (r(z)) is the displacement 
of the flux line at height z with respect to its average 
position in the {x, y) plane) 



(1.12) 



(ii) for T > T* , the potential well V(r) is shallow and 
U e (T) is strongly (exponentially) suppressed, 



U e (T)n^U (T/T 



)2 C -2(T/T*) 



(1.13) 



and the mean square width of the flux line is now very 
large 



{u\z)) ~ b 2 e 



2 2(T/T*) 2 



(1.14) 



A certain number of remarks seems to be appropriate 
at this uoint. First, it has been argued by Nelson and 
Vinokurl that in both regimes T < T* and T > T*, 
the flux line will be bound to the columnar pin, although 
only weakly so when T > T* . Altough this is perfectly 
true in the regime T <C T* where we expect many lo- 
calized states to exist and the physics to be dominated 
by these discrete, low lying states, at high temperatures 
T > T* , the problem is more delicate and requires care- 
ful attention. Indeed, when T > T* , we might fall into a 
situationllj where the only localized state is the ground 
state. In this case, the question as to whether the re- 
maining, extended states (belonging to the continuum 
spectrum of the quantum Hamiltonian TL) give a nonneg- 
ligible contribution to the partition sum (1.8) seems to 
be a perfectly legitimate one. If it turns out that the par- 
tition sum Z, in some physical regime, is dominated by 
those extended states, then we would expect this regime 
to correspond to a depinned flux line, i.e. a free line at 
high enough T 3> T* temperatures. 

Our second remark is related to the form of the mean 
square width of the flux line when T > T*. As T grows, 



it can be seen from equation (1.14) that (u 2 (z)) grows 



without bound and, assuming for the sake of argument 
an infinite superconducting critical temperature T c , at 
high enough temperature this last quantity migh become 
larger than the mean square width of a free flux line (here 
dj_ = d — 1 is the number of transverse dimensions) 



<« a (*)>o = 



d±k B TL 
12ei 



(1.15) 



This is, of course, conceptually unsatisfactory : the mean 
square width (u 2 (zj) of a "bound" flux line should be 



bounded from above by (it 2 (z))o, eq. (1.15), and, in the 
event that {u 2 (z)) is to become very large, its upper limit 
should not exceed (u 2 (z)} . This limit should not, how- 
ever, be put in by hand, but should emerge in a natural 
way from the formalism used to describe the physics. 

In what follows, we construct such a formalism. The 
key point will be to try to incorporate, in an ap pro ximate 
way, all the states {ip n } in the partition sum (1.8). This 
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is a rather difficult task within the Schrodinger formula- 
tion, since it requires the evaluation of a great number of 
eigenvalues E n of the quantum Hamiltonian TL and their 
associated eigenstates {ijj n (r)}. We thus have to look for 
alternate ways to solve our problem. 

For the purpose of achieving such an alternate, and 
hopefully more accurate, description of our system, let 
us rewrite the displacement vector r(z) as the sum 



y(z) = r + u(z) 



(1.16) 



where i"o is the position in the (xy) plane of the center of 
mass (CM) of the flux line, and u(z) is the displacement 
of the flux line at height z with respect to the CM posi- 
tion ro. The partition function Z of equation (1.3) can 
then be written in the form 



Jdr J Vu(z) exp { - /} J Q L dz [\ si (jf) * + V (r + u(z))] } 
I dr J Vv{z) exp { - f3 J Q L dz \h{^f) 



(1.17) 



where we have been careful to separate the CM from the 
internal modes, writing our measure of integration [dr(z)] 
as dr$ \du(z)]. Now let us imagine, just for the sake of 
argument, that we could succeed in integrating the above 
partition function over the internal modes {u(z)} exactW. 
The result would be a single integral over the CM modet^l 



Z 



I 



drr, 



-V cff (r )/k B T 



(1.18) 



where V e ff (r ) is the effective potential experienced by 
the CM mode after averaging over the internal degrees of 
freedom of the flux line. Since the original pinning poten- 
tial V(r) is very localized, we expect the effective poten- 
tial V e ff (r) to be localized as well. In fact, \V e ff(ro = 0)| 
would be nothing more than the exact value of the effec- 
tive pinning energy U e (T) per unit length of the flux line. 

An exact integration over the internal modes of the flux 
line being impossible to achieve, we have to resort to ap- 
proximations. An approximate but very accurate method 
which has been use d su ccessfully to calculate path inte- 
grals of the form (1.17) in various quantum statistical 
problems- is. the variational approach of Feynman and 
KleinertEJH. In this approach, one uses a Gaussian trial 
Hamiltonian to find an approximate expression for the ef- 
fective potential V e ff(r ) experienced by the CM mode. 
Applied |-tp various quantum problems (see for example 
Kleinerttj, and references therein), this approximation 
gives results which are in very good agreement with nu- 
merical simulations, both at low and high temperatures, 
which is a good indication that it does actually take into 
account, in a rath er a ccurate way, many more state s in 
the partition sum (1 



1. 



than the approximation fll.9| ). In 
the next section we apply this variational approach to 
our problem. In doing so, we shall only give the salient 
features of the calculation, and refer the,reader interested 
in technical details to Kleinert's bookll3, which gives the 
most complete and detailed presentation of the method 
present nowadays in the literature. 



II. VARIATIONAL APPROXIMATION FOR THE 
EFFECTIVE POTENTIAL 

We now introduce the following decomposition of the 
diplacement vector u(z) in Fourier modes q n = 2nir/L : 



u(z) = ^r( (Z „)e^ 2 



(2.1) 



The above decomposition is analogous to the decom- 
position of an imaginary-time trajectory in Matsubara 
modes, pihich is familiar from the equivalent quantum 
problemEd, or to the decomposition of a polymer's in- 
ternal mpdes into Rouse modes familiar from polymer 
physicso. Note that the n = mode is excluded from 
the summation since it corresponds to the center of mass 
mode ro. The Fourier coefficients r(q n ) are related to 
u(z) by 



r(<7» 



(2.2) 



In the ab ove r epresentation, the partition function Z of 
equation (1.17) can be written in the form : 



Z = I ^ /[,/u(:i] - / 







'dr\ 


/ dz 
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\dz) 



+ y(r + u(z))]} 



(2.3) 



where A is the transverse area of the system (in the plane 
which is perpendicular to the direction of the columnar 
pin). The integration measure [e?u(z)] stands forlia : 



[du(z)} = n 



dr re (q n ) dr lm (q n 



n=l 



hr/LPi-Lql 



(2.4) 



where r re (q n ) and Ti m (q n ) denote the real and imagi- 
nary parts of r(q n ), respectively, and where the denom- 
inator comes from integratin g the free path integral in 
the denominator, of equation ( 1.17 ) . Following Feynman 
and KleinertH-ffia, we introduce the following variational 
Hamiltonian : 



Hi 



1 



dz 



+n 2 (r )u 2 (z)] +£i(r„) (2.5) 
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where f2 2 (ro) and Li(ro) are unknown functions of the 
CM position ro to be determined by minimization of the 
variational free energy F v given by : 

F V = F 1 + (H-H 1 ) 1 

= F 1 + £dz [{V(r(z))) 1 - \ e 1 (^ 2 (r )u 2 (z)) 1 ] + 
- (LM) 1 (2-6) 

Here (--^i denotes an average taken with statistical 
weight exp(— Hi/ksT). On the other hand, F\ = 
— fcsTlnZi is th e free energy associated with the trial 
Hamiltonian (2.5), with 



) + 



f dr f. , , N1 r I r 1 r/dr 

Zi = j^Aj [du{z)]ex n-Tj dz 2 £i 

+ n 2 (r )(r(z)-r ) 2 ]}e- L ^/ T (2.7) 



We now use the Fourier decomposition (2.1) to rewrite 
the quantity /3-ffi in the form 

oo 

pHi = mi X>« + tt 2 (r )|r(<z„)| 2 + Li(r ) (2.8) 

71=1 

We thus obtain for the trial partition function Z\ : 

Zl =J ^. e -^(r ) J [du{z)] x 

OO 

x exp { - Lpi! Y,(<ll + fi2 ( r o)) [r 2 e (g„) + r 2 „ (<?„)]} 

71=1 

= J ^ e -^i(ro)$(n(r )) (2.9) 
Here and below, we denote by $(£l(r )) the quantity 



oo 

$(n(r„)) = n ( 



'17, 



XX \ q 2+Q2 {ro)/ 

m(r )/2 y-- 
sinh (m(r )/2) 



(2.10) 



x 2 \ sinhx 



where, in going from the first to the second-line, we used 
the well-known infinite product identitytZrta 

oo 



We now evalute the next term in the trial free energy 
F v . We have : 

(V(r(z))>i = Zf 1 / ^ e-? 1 *™ [d[u(z)]e-f>* x 



V(k) e 



k -[ ro+ E„vo r ('?") el '""] 



(2.11) 



where the Fourier transform V(k) is given by 



V(k) = J dr ^(r)e" ik 



(2.12) 



and where we used the shorthand notation J k = 
/ d d± ]s./ (2ir) d± . Writing the sum J2 n ^o r (ln) m the form 
(here and below, c.c. denotes complex conjugation) : 

oo 

]T r(q n )e^ z = £ (r(9„)e l9 " z + c.c.) 



n#0 



2 [iY e (<7„) cos(q n z) - r im (q n ) sm(q n z 



and perfo rmin g the resulting Gaussian integrations in 
equation (2.11) above, we obtain : 



(V(r(z))) 1 = ±J ^e-^M$(Q(r ))F(r ) (2.13) 



In the above expression, we defined the function 

-ifc 2 a 2 (r )+ik-r 



V(r ) = / V(k) e"*' 



(2.14) 



which can be interpretectJ as an average of the origi- 
nal pinning potential V(ro) over a region of radius a(i"o) 
around r , where a(r ) is given by : 

2 °° 1 

fl2(ro) = Z^g^TT^R 

cotn 1 ' 



L/3iiil 2 (r )\ 2 



where, in going frqia the first to the second line, we made 
use of the formula 



El 7T 1 

— 5 o = 7T~ coth(7ra) - — 
n 2 + a 2 2a 2a 



(2.16) 



The third average in the trial free energy fl2.6| ) is eval- 
uated in a similar fashion. We have : 



n 2 (r )u 2 (z)) 1 =2j2(n 2 (r )\r(q n )\ 2 } 1 



1 f dr Q 



fi 2 (r )e- /3il(ro) / [du(z)] x 



ZiJ A 

OG 

exp | - L^ex ^(q 2 + n 2 (r )) [r re (qi) + r irn ( qi )] j 



2 5Z ( r re(9n) + r L(9")) 



(2.17) 



Again, evaluating the Gaussian integrals over the internal 
modes r(q n ), we obtain : 
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dz -£! (tf(r )u 2 (z)) 1 



ZiJ A 



e -/3Li(r ) x 



x -d ± Ze x $(fi(ro))n 2 (r )a 2 (r ) 



(2.18) 



The last average we need for the evaluation of F v is 
(Li(i"o))i, which can be easily obtained, and is given by : 

(ii(ro)>i = yJ^ ^(^(ro))e- /3il ^Li(r ) (2.19) 

Collecting all terms, we finally obtain for the trial free 
energy F v the following expression : 



F v =F 1 + ±-J ^° *(n(r )) e 



•/8ii(ro) 



x [LF(r ) - -djLEifi 2 (r )a 2 (ro) - L x (r )] (2.20) 

Variation of the above free energy with respect to LjI 
shows that F v is extremal, and is in fact minimalo 
for 

Li(t ) = LV(r ) - ldj_Leifi 2 (r )a 3 (ro) (2.21) 
The variational free energy F v thus reduces to : 

F V =F 1 = -k B T\nZ 1 (2.22) 



with Z\ given by e xpression (2.9), with Li(ro) replaced 
by the result ( [2.21 ) above : 



z ^ = I d£0 $(Q( ro )) e -/9[iV'(>-o)-^ii£i!! 2 (ro)a 2 (ro)] 

■A. 



From this last equation, we see that Z v = e 



-0F V _ „-/3Fi 



can be written in the form 

Z», - 



dr o p -0W(r o ) 



A 



(2.23) 



This is exactly the goal that we set out for in the intro- 
duction : W(ro) is no more than the effective potential 
experienced by the center of mass of the flux line when 
all the internal modes have been integrated out. Within 
our variational approximation, it is given by : 

W(r ) = -k B T In <P(n(r a)) +LV(r ) + 

- ^ x £ei0 2 (r )a 2 (ro) (2.24) 

Furthcx jiainimization with respect to f2 2 (ro) leads to the 
resultlli'EJ that this last quantity is given by 



dV(r ) 



dx£i da 2 (r ) 

-1- / k 2 V(k)e~ 
d±ei J k 



jfc 2 a 2 (r )+ik-r 



(2.25) 



Taking for the pinning potential of a columnar pin the 
Gaussian form 



V(r) = -U exp(- 



2b 2 ' 



2(6 2 + a 2 (r )) 



we find for V'(ro) the following expression 

2 

V(r ) = -Ux exp ( 
where U\ is given by : 
Ut = U ( 



68 



b 2 + a 2 (r ) 



d ± /2 



(2.26) 



(2.27) 



(2.28) 



In the following we shall mostly be interested in the ef- 
fective depth of the pinning potential, and hence we shall 
be interested only in the values of il 2 (r ) and a 2 (r ) eval- 
uated at ro = 0. Using th e exp ression ( 2.26 ) of the pin- 
ning potential in equation ( 2.25| ) and performing the k in- 
tegral, we easily obtain the following result for f2(ro = 0) 
(henceforth, we shall take d± = 2) : 



n(o) 



ei 6 2 + a 2 (0) 



(2.29) 



At very large sample thickness L and low enough tem- 
perature, the l owe st lying, localized states dominate the 
partition sum (1.8). We expect in this regime a localized 



flux line, with finite values of the parameters J7(0) and 
a(0). Assuming that Q(0) remains finite (nonzero), so 
that Ltt(0) is very large when L — ► oo, we can approx- 
imate a 2 (0) from equation ( p. 15 ) in the above limit by 
the following expression : 



« 2 (0) 



1 



2/3ei^(0) 



DC 



(2.30) 



Using equation (2.29) to eliminate O(0) from equation 
fl2.30| ), we obtain the following equation for a 2 (0) : 



a 2 (0) 



b k B T 



a 2 (0) 



b 2 o 



(2.31) 



Solving the above equation for a 2 (0), we obtain the re- 
sult : 



a 2 (0) 



boksT 



1 



k B T 



(2.32) 



which shows that the mean square width of the flux line 
a 2 (0) diverges at the temperature Td p given by : 



ksTdp = 2b J Uq ii(T dp ) 



(2.33) 



Before proc eedin g any further, we need to observe that 
expression ( 2.32Q for a 2 (0) is valid only in the tempera- 
ture range < T < Td p (for T > Td p , the denominator 
in equation (2.32) becomes negative, which would violate 
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the requirem ent a 2 (0) > 0). We therefore expect the ap- 
proximation ( [2.30 ), which relied on the assumption that 
f2(0) was finite (i.e. nonzero), to be no longer valid when 
T > Tdp. Assuming, in this regime of temperatures, that 
Lf2(0) <C 1, and using the Taylor expansion of the hy- 
perbolic cotangent near x = : 



1 x 



cothx = — I h o(x ) 

x 3 



x -> 



we obtain 



m , /LCi\ L 2 n 2 (o) 

coth - 1 ~ — 

2 V 2 / 12 



Inserting this expression back into equation ( 2.15| ), we 
obtain : 



/ 2\ j 2/n^ d±k B TL 
(u ) = d ± a (0) = — 



(2.34) 



which is nothing more t han the mean square width of a 
free flux line, equation ( 1.15 ). 

The above results suggest that the flux line undergoes 
at the critical temperature T = T c i p a transition from 
a low temperature (T < Td p ), pinned phase to a high 
temperature phase (T > Td p ) where the flux line is delo- 
calized. Note the si milari ty of our "depinning" temper- 
ature Tdp, equation (2.33), to the^-characterist ic tem per- 



ature T* of Nelson and Vinokumu, equation ( 1.11 ); the 



two expressions differ only by a constant of order unity, 
which should in fact be considered good agreement, given 
the fact that different analytic forms for the pinning po- 
tential are used in the two cases. Unlike the quantum 
approach however, our variational met hod y ields a mean 
square width (u 2 ) = d±a 2 (0), equation (2.34), which does 



not gr ow without bound with temperature as in equation 
(1.14). The fact that we are able, in the high tempera- 
ture phase, to obtain the correct expression for the mean 
square width of a free line from our equations is an in- 
dication that o ur va riational approach is superior to the 
approximation (1.9), and that the resulting scenario of a 
true phase transition might indeed be the correct way to 
describe the physic s of ou r syste m. 

Equations &2$ ), ([Eg|)-(|!2§), and (p|)-(p|), are 
t he m ain results of this paper. In particular, equation 
( [2.24 ) gives the effective potential experienced by the cen- 
ter of mass of the flux line after all the internal modes 
have been integrated out, and as such, it allows us to dis- 
cuss the physics of the flux line near the columnar pin in 
terms of only one variable, namely the position ro of the 
center of mass mode. This discussion will be the subject 
of the next section. 



III. DISCUSSION AND CONCLUSIONS 

We now can construct a complete and consistent pic- 
ture of the physical behaviour of a single flux line near 



a columnar pin. As we did in the introduction, we shall 
consider two separate temperature regimes : 

(i) For T < Tdp, the internal fluctuations of the flux 
line smear out the pinning potential only weakly, so that 
the center of mass of the flux line still experiences a fi- 
nite (nonvanishing) effective potential W(tq), We can 
therefore think of the center of mass mode of the flux 
line r as an ordinary classical particle trapped in the 
pinning potential W(tq), and in contact with a ther- 



mal reservoir at temperature T. Using equations (2.24), 
( j2~27| ), ( gTgg ) and ( pT32[ ), and the definition Q2.10|) of 
the function one can easily show that the strength 
U t (T) = \W(r = 0)| of the effective potential W(r ) 
(which is also the total pinning energy of the flux line) 
can be written in the form Ut{T) — LU e (T), where the 
effective energy per unit length U e (T) 



U e (T) = U (l 



k B T 



2W C/ o£i(T) 



(3.1) 



is reduced with respect to the bare pinning energy per 
unit length Uq and \fajpishes exactly at T = Tdp. Since 
the average exit timai§tj r e of a classical (overdamped) 
particle in contact with a thermal bath at temperature 
T, from a potential well of depth Ut is proportional to 
the Boltzmann factor exp(Ut/ksT), we see that, in the 
case of the center of mass mode of the flux line where 
exp(U t {T)/k B T) = exp(LU e (T)/k B T), no matter how 
small U e (T) is, in the limit of a large sample thickness 
L, r e will be very large, so that ergodicity is effectively 
broken and the flux line may be considered, for all practi- 
cal purposes, as localized in the vicinity of the columnar 
pin. In the boson language, this is what is meant by say- 
ing that the flux line is "bound" to the columnar defect. 
In this low temperature regime, only low lying, localized 
states give a si gnifi cant contribution to the sum flLS| ), 
and expression (L9) is "good enough" an approximation 
to the partition function. 

(ii) For T > Td p , the internal modes fluctuations di- 
verge, and the mean square width of the flux line takes 
on its "free" value. As a result, the pinning potential 
is completely smeared out : the effective pinning energy 
Ut(T) = | W(ro = 0)| experienced by the center of mass 
mode v anish es (as c an eas ily b e verified using equations 
(|2 lOP , Q2.24D , ( |2 29D and fl2.34|)), and the flux line is not 
pinned at all by the columnar defect. In this regime, 
(if it does actually exist, and is not just an artifact of 
our variational approximation) , one has to be careful not 
to neglect the extended states (in the continuum part 
of the spectru m of the quantum Hamiltonian TL) in the 
partition sum (1.8), as these states turn out to dominate 
the high temperature behaviour. The flux line wanders 
freely around the columnar pin, and this free behaviour 
emerges in the most natural of ways from our equations. 

In conclusion, in this paper we applied the powerful 
and yet simple Feynman-Klcincrt variational method to 
investigate the behaviour of a flux line near a columnar 
pin in type II superconductors. Our results suggest that 



G 



the flux line undergoes a phase transition from a low tem- 
perature phase where it is effectively pinned to the defect, 
to a high temperature phase where the line behaves as a 
free flux line. While there are some good indications t\ 
our theory might be more accurate than previous wort 
the issue of whether a true phase transition does indeed 
exist remains open. It therefore appears that the only 
way to decide whether there is an actual phase transi- 
tion at T — Td p is to resort to other techniques such as 
renormalization group methods or numerical simulations. 
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IV. APPENDIX A : FLUCTUATIONS OF A 
SINGLE FREE FLUX LINE 

In order to make this paper self-contained, we here 
present the derivation of the mean square wid th (it 2 ) of 
a single free flux line. The Hamiltonian (1.1) of a free 
flux l ine can be rewritten within the Rouse decomposi- 
tion (2.1) in the form 



H = £ ULix) q 2 n \r(q n )f 



(4.1) 



From this expression, we see that the free elastic propa- 
gator Go(q n ) such that 



(r a (q n )rp(q m )) = k B TG (q n )S a ,f]O~n,- 
is given by 

G {q n ) 



(4.2) 



(4.3) 



Le x ql 

so that the mean square width of the flux line is given by 



(u 2 (z)) =d ± k B Tj2Go(qr, 
d±k B TL 



(4.4) 



where, in going from the first to the second line, we used 



the fact that J2 r . 



7T 2 /6. 



For completeness, we here also compute the average 
"tipping angle" ((dr/dz) 2 ) of the free flux line as well as 
the relative displacement ([11(2) — u(z')] 2 ). The average 
tipping angle of the line is given by 



((dr/dz) 2 ) = 2d 1 _k B TY j qlGo{qn) 



(4.5) 



where, in order to obtain a finite result, we introduced 
the cut-off ~ L/£o which corresponds to summing 
over all wave- vectors q n such that \q n \ < 2ir/£ c (the fac- 
tor 2 on the rhs of the above equ ation comes from the 
negative q z modes). Equation (4.5) then yields 



((dr/dz) 



2d±k B T 



(4.6) 



which is -basically the result quoted in Appendix B of 
referencetl. In the regime of temperatures where the co- 
herence length £, C (T) along the easy axis of a layered 
material is smaller than the layer spacing s, we cut off 
the sum in equation ( f4.5j ) at \q n \ < 2ir]s, which gives 
((dr/dz) 2 ) = (2d±k B T/i!s) instead of (Q. 

For the mean square displacement ([u(z) — u(0)] 2 ) we 
find : 



<[u(z) - u(0)] 2 ) = 4d x k B Tj2 G o(<?n) [1 - cosq n (z)} 



Using expression (4.3) of the elastic propagator, and 
transforming the sum into an integral, we obtain 



(K*)-u(o)] 



2d±k B T f 00 1-cosgz 
" ' ' dq 



11 Si 



q 



Changing the variable of integration form q to u = \z\q, 
we finally obtain 



([u(z)-u(0)] 2 ) = 



d I fcnT 



where we used the fact thaiS J °° du (1— cos u) /u 2 



(4.7) 
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